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We examine the effect of including the 'combined' HERA structure function data in the MSTW 
global fit for parton distribution functions (PDFs). The combined neutral-current HERA data have 
a significant, if not dramatic, effect, of up to 2-3% at NLO for Z boson and Higgs production at 
the Tevatron and LHC, and a generally slightly smaller effect, particularly on LHC processes, 
at NNLO. This is an amount comparable, or less than, the typical PDF uncertainties, and hence 
we do not intend to release an imminent update to the MSTW 2008 fit. We also investigate the 
consistency with the recent D0 data on electron and muon charge asymmetry from W decays and 
the direct CDF measurement of the W charge asymmetry. The D0 lepton charge asymmetry data 
imply a fairly large change to the down-quark distribution and/or large nuclear corrections to be 
applied when fitting to deuterium structure function data, while the CDF W charge asymmetry 
data are more consistent with the existing PDFs. However, it is difficult to reconcile all of the 
Tevatron W —?' iv charge asymmetry data sets with the fit, and to some extent, with each other. 
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1. HERA combined data 

The MSTW 2008 global fit ^ for PDFs contains a very large variety of data sets, including 
a number from HI and ZEUS on structure functions. These structure function data (along with 
some newer sets) have recently been combined in [Q]. The increase in precision comes about not 
only from a combination of statistics, but from the fact that one collaboration often controls a 
source of systematic uncertainty better than the other, so the systematic uncertainties can be greatly 
reduced by combination. In addition, this treatment of the correlated errors means that the central 
values are not simply the weighted averages. These data were fit in ^ and differences between 
previous HI and ZEUS fits noted. Here we, too, investigate the inclusion of the combined neutral- 
current (NC) data instead of the component sets, adding the errors in quadrature for the moment 
(systematics in the combined data set are small). We include the uncombined charged-current data 
which are statistics dominated. We fit to data with >2 GeV^ (553 pts.), but also look at results 
for 2^ > 3.5 GeV^ (524 pts.) to compare with All other details of the fit are as in |jl|]. 

At NLO the fit quality achieved is 
2610/2471 for the total data. For the HERA 
NC data it is 600/553 and 530/524, com- 
pared to 483/524 in [Q]. In order to check 
if the worse quality is due to the other 
data in the MSTW fit we also fit only to 
HERA structure function data. This results 
in 515/553 and 445/524, now much bet- 
ter than the HERA fit [^, presumably due 
to extra parameterisation freedom. Clearly 
the fit quality is significantly affected by 
tension with other data sets. The MSTW 
global fit with HERA combined data requires 
asiM^) = 0.1215, a little higher than the de- 
fault 0.1202. Keeping at the default value re- 
sults in a fit quality 10 worse for HERA data 
and the global fit. Our 'HERA data only' fit 
gives as{M^) = 0.123. The resulting PDFs 
are shown in Fig. [I] as a ratio to MSTW08 
with the la uncertainty. For the global fits 
the up quark strays outside this uncertainty 
at X ~ 0.01, otherwise there is little change. 
The predictions for hadron collider processes give 2-3% variations for Z production, but less for 
Higgs production (from gluon-gluon fusion, with Mh = 120 GeV), see Table [|. Our fit to 'HERA 
data only' gives a similar small-x gluon but for x > 0.01 it is far softer, while the up-quark dis- 
tribution is not well constrained in shape for x > 0.01. These PDFs compare very badly to many 
unfitted data sets. Both the variation and comparison to unfitted data are far worse than for the 
PDFs in showing that it is implicit constraints on the quark parameterisation rather than a real 
data constraint that render these PDFs similar to those from global fits at high x. 




Figure 1: The ratio of NLO PDFs fit to the combined 
HERA data to the default MSTW08 distributions. 
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At NNLO the fit 
quality for the HERA 
NC data is 575/553 and 
529/524, better than NLO. 
0:5 (M|) moves only from 
0.1171 to 0.1181, and 
fixing it leads to ~ 
8 higher. Fitting only 
HERA data gives 494 /553 
and 467/524. The change 
in PDFs is similar to 
NLO, but with a tendency 
to dip slightly at;c < 0.001. 



PDF set 


fi/+/--Oz(nb) anipb) 


B/+/--az(nb) Of/Cpb) 


NLO 


Tevatron 


LHC (14 TeV) 


MSTW08 
New HERA 

fix as (Ml) 


0.243+lfZ" 0.746+^°^" 

-2.0% -4.4% 

+3.1% -0.7% 
+3.0% -4.0% 


2.00+?,lJ 40.7+^"J 

^"-1.8% -2.7% 

+2.5% +1.2% 
+1.8% -0.8% 


NNLO 


Tevatron 


LHC (14 TeV) 


MSTW08 
New HERA 

fix a5(M|) 


n 0*51+2.2% QQj;j;+5.4% 
"•^-'^-1.8% "•^-'-'-4.8% 

+3.0% +0.2% 

+2.6% -2.9% 


o n'^+2.6% cf^ cr+3.6% 

z.UJ_2.1% -'"•-'-2.7% 

+1.8% +0.8% 
+1.2% -1.0% 



Table 1: The cross sections at the Tevatron and LHC with "PDF+as" un- 
certainties and the changes in fits using the new combined HERA data. 



The NNLO predictions have generally less variation at the LHC than at 
NLO, see Table [l[ The effect of the combined data is significant, and will be included in an updated 
set soon, but certainly is not dramatic enough to invalidate the present MSTW 2008 PDFs ||l|]. 

2. Tevatron lepton charge asymmetry from W decays 

There are new D0 data on electron [Q] and muon [Q] charge asymmetry and CDF W charge 
asymmetry data [||]. These are far more precise than the previous measurements [§, 0], which 
already give the main constraint on some PDF eigenvectors. These new data should constrain the 
down quark for 0.01 <x < 0.7, where the current main constraint is deuterium DIS and is subject 
to uncertainty from nuclear corrections (a source not included in PDF uncertainties). The current 
fit to asymmetry data is moderate At NLO = 25/10 (D0) [|] and x^ = 29/22 (CDF) [|]. 

The new data cause 
worse problems. Stan- 
dard MSTW fits give a 
very poor comparison to 
both D0 electron and 
muon data, as seen in Ta- 
ble ^ and Fig. ^. We have 
tried a variety of alterna- 
tives, weighting the asym- 
metry data and/or making 
cuts on other data in the 
fit. The fit quahty to D0 e 
and jU data and other data 
for these variations is also 
shown. Here 'cut' means 
the omission of BCDMS 
proton and deuterium data 
and NMC n/p data which 
are very badly described. 
NNLO corrections | 



fit 


zVl2(e) 






ZV2689 


zVl6(At) 


PT (GeV) 


>25 


25-35 


35-45 


non- D0 


two bins 


Standard 












2008 data 


116 


19 


144 


2518 


542 


D0, 


71 


23 


81 


2551 


358 


D0,(w) 


25 


10 


23 


2942 


183 


D0^ 


26 


55 


58 


2640 


119 


D0^ (w) 


33 


79 


88 


3131 


10 


D0^ cut 


33 


52 


55 


3190 


26 


Deut. Corr. 












2008 data 


25 


32 


42 


2455 


140 


D0,(w) 


25 


9 


23 


2551 


192 


D0^ (w) 


38 


67 


75 


2649 


11 


^0e+ti (W) 


24 


16 


40 


2848 


42 


D0^ Pt>25 


23 


38 


32 


2454 


229 



Table 2: Description of D0 lepton asymmetry and remaining data, without 
(upper) and with (lower) deuterium corrections, (w) denotes a high weight. 
Asymmetry data included in each fit are shown in bold type. 



, 10] do help the fit, but only marginally. 
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0.3 



> 25 GeV, > 25 GeV 



0.2 



0.1 



S -0.1 



-0.2 



-0.3 



-0.4 



-0.5, 





MSTW08: = 278 (8 pts.), = 116 (12 pts.) 
Fit new D0 '\: = 30, x' ='26 
Weight by 100: x^ = 5, x'= 33 
Cut BCDMS4-Ni«Cn/p: y}'= 4, = 33 
Deut. corr., fit oid D0 A^f: x^ = 48, x^ = 25 
Deut. corr., fit new D0 A, 
Deut. corr., fit new D0 A^: x' 
D0 (prei.) A , L = 4.9 fb ' 
D0 (pubi.) a" , L = 0.75 fb 



0.5 



I . 
1.5 



2.5 



> 35 GeV, > 25 GeV 



In order to try to improve the 
fit, an extra parameter was added to 
both the high-x valence quark dis- 
tributions. This had negligible im- 
pact in the fit quality and extracted 
PDFs. Then variations in the deu- 
terium corrections to the structure 
function data fit were tried. In the 
standard MSTW fit these data have 
small corrections for shadowing at 
small-x, but none at high-x. Re- 
moving these corrections and refit- 
ting using the standard MS TWOS 
data leads to x'^ = 19/10 (D0) 
and = 25/22 (CDF), a signif- 
icant improvement. The up and 
down quarks change by 1-2% near 
X = 0.02. Given this mild success 
we also tried a more sophisticated 
approach to corrections for deu- 
terium data, i.e. 2^ -independent 
deuterium corrections for all x ap- 
plied to theory by means of a 
smooth function with 4 free pa- 
rameters. This improves the qual- 
ity of the fit to non-asymmetry data 
significantly, as seen in the lower Figure 2: The fit to D0 lepton asymmetry data for all pj > 
half of Table |, and using the stan- 25 GeV (top) and 35 GeV <pt< 45 GeV (bottom), 
dard MSTW08 data sets gives x^ = 6/lO (D0) and = 21 /22 (CDF). When fitting to the newer 
D0 asymmetry data the free deuterium corrections also help a great deal, as seen in Table |2| and 
Fig. |2[ A good fit (given the scatter of data) can be found for the electron data in the combined pj 
bin, with no deterioration in fit quality for other data types. However, although very good fits to the 
electron data in separate-/??- bins, or to the muon data, can be found, they both, especially the latter, 
result in a deterioration in the fit quality to other data. We also find that the muon data and elec- 
tron data cannot be fit near their best individual quality simultaneously, see Table ^ The required 
deuterium corrections are shown in Fig. |3[ The general shape is similar to expectations, labelled 
"simple model", but in all cases, especially when fitting muon data, the correction is low in the 
region of x = 0.1. The main change in the PDFs, when deuterium corrections are applied, is an in- 
crease in the d{x, Q^) for x > 0.02, and, in particular, a 10% increase for x 0.4 at = 10^ GeV^. 

As well as tensions between the two D0 data sets there are also conflicts with the quality of the 
comparison to CDF W asymmetry data . The MSTW08 PDFs give (an approximate) « 28/ 1 3 
which, given the scatter of points, is quite good. The MSTW fit to the standard data with deuterium 
corrections gives x^ ~ 24/13, and as seen this also fits the combined-/??- D0 electron data well. 
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Good fits to the separate pr electron data and/or muon data give > 100 for comparison to the 
CDF W asymmetry data, due to too much asymmetry. i < 

We conclude from this study that inclusion of deu- 
terium corrections and an investigation of their uncer- 
tainty is important for global fits. The deuterium data 
can be fit without them, but they improve the com- 
parison even with the older low statistics asymmetry 
data. An examination of the more recent asymmetry ^ 
data leads us to conclude that the maximally consis- 
tent sets in the fit are the CDF W asymmetry data and | 
the combined-7?7- D0 electron data. However, a good g 
fit to the latter, without seriously affecting the rest of 
the global fit, requires slightly unlikely deuterium cor- 
rections. In summary, it seems at present that the diffi- 
culties in reconciling the different asymmetry data sets 
with theory, and with each other, necessitates further 
study, both of data and the theory to be applied, before 
the true impact in a global fit can be understood. 




Simple model 
Low liimi data 
DOII muon 
DOII electron 
DOII electron combined 
DOII electron combined weighted 



Figure 3: The deuterium corrections. 
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